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INTRODUCTION 

FIESTA ROC (Finite Element Semiconductor Three-dimensional Analyzer by Ralph 0. Clark) is a computational tool 
for investigating in detail the performance of arbitrary solar cell structures. As its name indicates, it uses the finite element 
technique to solve the fundamental semiconductor equations in the cell. It may be used for predicting the performance 
(thereby dictating the design parameters) of a proposed cell or for investigating the limiting factors in an established design. 

THE FINITE ELEMENT METHOD 


The fundamental semiconductor equations solved by FIESTA ROC are [1] 


V . (e£) — p (Poisson’s Equation) 

(la) 

V ♦ T n = G — R (Electron Continuity Equation) 

(lb) 

V • T p = G — R (Hole Continuity Equation) 
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Here £ is the electric field, e the permittivity, p the volume charge density, G — R the net carrier 
minus recombination), and T n and T p the electron and hole flux densities. In turn, £, T a, 
electrostatic potential ip and the quasi-Fermi potentials <p n , (p v by 
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T p — V(f> p 

(2c). 


As boundary conditions, we assume that the domain of simulation, fi, has a boundary dQ that can be partitioned into 
two segments: the Dirichlet boundary on which ip, <p n , and <pp are prescribed, and the Neumann boundary <9£2jv. on 

which the components of V^\ V0 n , and V</>p normal to the surface all vanish. In particular, on the Dirichlet boundary, 
which corresponds to the contacts of a device, an infinite surface recombination velocity is assumed, pinning n and p to their 
levels at thermal equilibrium and <p n and <p p both to ip — ipo. On the Neumann boundary, arbitrary surface recombination 
velocities are possible. Note that, for a one-dimensional model of a two-terminal device under these assumptions, we cannot 
model finite surface recombination velocities at the ends of the device because both ends represent contacts, i.e. Dirichlet 
boundaries. 

Equations (1) can all be written in the generic form 

Vf - s = 0 (3). 

To solve (3) by the finite element method, we first write it in its weak form, as follows. Let w be an arbitrary piecewise- 
continuous function on Q, vanishing identically on d£lo Multiplying (3) by w, integrating over the entire volume, and 
integrating the first term by parts, we obtain 

0 = / w(V • f - s) = I rt f n - f (Vw f + ws), 
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where u is the unit outward normal to the surface. Note that since w vanishes on the Dirichlet boundary and f ■ n on the 
Neumann boundary, the surface integral vanishes, leaving us with the weak form of equation (3): 



u; ■ f + w s) ~ 0 


( 4 ). 


We now discretize the problem by fixing N points, called nodes, between which w, 4> n , and d>„ all vary in a 
n P art ' cu * ar > we express <f>„, <p p ] as a linear combination of continuous, piecewise smooth functions 
tv.fY'.VP, ,<p i J, I _ 1, N, each of which assumes the value of 1 at node t and zero at all other nodes Thus 

" = 1 W fr t= J-* = EL *?*?. 4> P = *?. where the are arbitrary coefficients (to keep w 

arbitrary) and V’.'. 4> % . </>, are just the values assumed by ip, <j> n , and <p p at node i. Now equation (4) becomes 


J2 Wi / (VW<-f+ WJ*) = 0 

.=i 

Since the N coefficients w { are arbitrary, equation (5) really represents the system of TV equations 

/(VW, f + Wis) = 0, i = N 
Jn 


(5). 


( 6 ). 


In fact with our three coupled equations (1) and three dependent variables * <j> n , 0„, we now have 3 N coupled non- 

Imear equations ,n the 37V unknowns , <f ,\ , i = 1 , . . . , AT. These are solved by the finite element code by means of a 

generalized Newton s method. 


ADVANTAGES OF FINITE ELEMENT 

B ecause the finite element method is an integral method, with only first derivatives appearing, many complications 
associated wth boundary value problems disappear. First of all, surface efFects. whether at a boundary between regions 
o he device (e g the junction) or at the physical surface of the device (e.g. surface recombination) are effectively a 
special case of bulk efFects where the integrand includes a Dirac delta function to localize the effect to the boundary The 
delta function converts volume integrals into surface integrals. Thus, surface recombination and interface charge density 
r ? a ^. .V" * eaSI y y surface integrals in equation (6). In addition, the perennial problem of matching the 

electric field at a jump discontinuity of c does not even come up, having been integrated out of existence. Thirdly, multiple 
non-interacting recombination levels are dealt with simply by including additional terms in R (which appears in the s in 
equation (6)). Finally, non-rectilinear elements are readily handled by performing the integrations in (6) numerically 


IMPLEMENTATION 

r W * have ^P ,emented a one-dimensional prototype of FIESTA ROC in C, running on an IBM-PC compatible and on a 
Lray XMH. I he prototype features automatic mesh generation and automatic dark and light I-V simulations, as well as 
spectral response. Arbitrary numbers of recombination levels, degenerate statistics, and heavy doping effects are supported 
Finite surface recombination velocities, as noted earlier, are imcompatible with our choice of boundary conditions in the 
one-dimensiona case, so they have been left out of the prototype. Once the solution to equation (6) has been found for a 
given b'as and illumination, many quantities of interest can be examined at any point in the cell. These include the electron 
and hole current densities, electric field, electric potential and quasi-Fermi potentials, net recombination, n, and p. 


EXAMPLE 

As an example of the use of FIESTA ROC, we have modeled an InP shallow homojunction n+pp+ cell from Spire 
Corporation which has been extensively studied with analytical models [2], The growth parameters and modeling parameters 
ound in [2] are reproduced in Tables I and II. Without any change in the parameters of either table, the results of light I-V 
and spectral response simulations by the prototype are shown in Figures 1 and 2. In addition, the measured and calculated 
solar cell parameters are shown in Table III. 
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Tabic I: Growth Parameters for Spire 6 Cell 


emitter width 

400 A 

n+ emitter doping (Si) 

1 x 10 18 cm” 3 

base width 

3 pm 

p base doping (Zn) 

2 x 10 18 cm" 3 

BSF width 

250 pm 

p + BSF doping (Zn) 

5 x 10 18 cm" 3 

cell area 

0.25 cm 2 

grid shadowing 

4.8% 


Table II: Modeling Parameters for Spire 6 Cell [2] 


indirect lifetime of holes in emitter 

2.0 ns 

hole mobility in emitter 

75 cm 2 /V-s 

indirect lifetime of electrons in base 

150 ns 

electron mobility in base 

3988 cm 2 /V-s 

indirect lifetime of electrons in BSF 

0.60 ns 

electron mobility in BSF 

2456 cm 2 /V-s 


The real utility of a numerical solver, however, lies in its ability to show what is happening in various regions of the 
cell. In this respect, it functions as a sort of computational microscope, allowing the investigator to probe all regions of 
the cell to determine, for example, which contribute most to the total recombination, light-generated current, and so on. 
In Figure 3, we show the calculated net recombination for three different bias points — short circuit, maximum power, and 
open circuit — plotted against the spatial coordinate x. Here x = 0 represents the metallurgical p-n junction and x = 3/xm 
the base-BSF low-high junction. The emitter, being only .04/im wide, is invisible on this scale. The junction space charge 
region is visible as a low- recombination valley at short circuit and a high-recombination peak under forward bias. Note, 
however, the additional peak in recombination in the BSF region. In Figure 4, we show the same curve, at maximum power 
only, calculated for the same cell and for another differing only in the doping in the BSF, namely 5 x 10 17 rather than 
5 x 10 18 cm -3 . The indirect recombination lifetime in the BSF is also correspondingly higher (6.0 ns instead of .60). The 
peak in the BSF almost disappears when the BSF doping is decreased. Such an analysis can aid cell fabricators in designing 
and producing more efficient cells. 


CONCLUSIONS 

We have demonstrated a one-dimensional prototype of a flexible finite element package that will enable cell designers to 
simulate a variety of effects and to pinpoint problems in proposed or existing cells. The two- and three-dimensional versions 
are currently under development. 
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Figure 2: Spectral response for cell Spire 6, with calculated results by FIESTA ROC prototype. 
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Figure 3: Net recombination calculated by FIESTA ROC prototype at three bias conditions. The n+p junction is at c 
and the pp + junction at x = 3pm. 
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Figure 4: Net recombination calculated at 
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iax power voltage for cells with different BSF dopings. 


31-6 



